Quantum Adiabatic Computation With a Constant Gap is Not Useful in One 

Dimension 
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We show that it is possible to use a classical computer to efficiently simulate the adiabatic evolution 
of a quantum system in one dimension with a constant spectral gap, starting the adiabatic evolution 
from a known initial product state. The proof relies on a recently proven area law for such systems, 
implying the existence of a good matrix product representation of the ground state, combined 
with an appropriate algorithm to update the matrix product state as the Hamiltonian is changed. 
This implies that adiabatic evolution with such Hamiltonians is not useful for universal quantum 
computation. Therefore, adiabatic algorithms which are useful for universal quantum computation 
either require a spectral gap tending to zero or need to be implemented in more than one dimension 
(we leave open the question of the computational power of adiabatic simulation with a constant gap 
in more than one dimension). 



There arc many different models for quantum compu- 
tation. The most standard approach is the gate model, 
combined with appropriate error correction to deal with 
decoherence^l]. Other approaches include measure- 
ment based quantum computer 0, topological quantum 
computing [ij], and adiabatic quantum computation 
Adiabatic quantum computation is very natural because 
one can imagine slowly changing the Hamiltonian follow- 
ing a path in parameter space, starting from some simple 
Hamiltonian with a known ground state, and arriving at 
some final Hamiltonian whose ground state encodes the 
solution of a difficult optimization problem^. 

While adiabatic quantum computation has been shown 
to allow for universal quantum computation Q, and hence 
is equivalent in its computational power, the fault toler- 
ant approach and related threshold theorems Q have not 
been generalized to adiabatic quantum computation. In- 
stead, one can rely on the spectral gap in the Hamiltonian 
to protect against errors. The spectral gap is also in- 
teresting in adiabatic quantum computation because the 
time required to perform the computational scales with 
the inverse spectral gap. It has at least been shown [7J 
that one can produce a constant gap against local noise 
on 1 and 2 qubits, but has never been shown that one 
can produce a constant gap against all excitations. 

Therefore, it is of great interest to determine if univer- 
sal adiabatic quantum computation can be performed in 
systems with a spectral gap of order unity. In a sense, 
topological quantum computation provides a means of 
performing universal adiabatic quantum computation 
with constant gap, by adiabatically changing the Hamil- 
tonian to drag defects around each other. However, this 
topological approach relies on having a large ground state 
degeneracy. In this paper we show that, at least in one 
dimension, adiabatic quantum computation in systems 
with a unique ground state and a constant spectral gap 
is not useful for quantum computation as it can be sim- 
ulated efficiently on a classical computer. 

Main Result — We consider the following problem. 



Consider a parameter-dependent Hamiltonian H(s) — 
5Zj hi,i+i(s), with hi t i-\-i(s) having support on sites i, i+1 
and with < J, so that interactions are near- 

est neighbor. Let there be N sites, and assume that 
each site has a Hilbert space dimension D which is 0(1). 
Assume that for all s with < s < s max we have a 
spectral gap AE with J/AE being 0(1). Finally, as- 
sume that ||9 s ft,j ; j_|_i|| < J; this last requirement simply 
sets some scale for how large s is. In general, if we have 
H^s/i^i+ill < X, for any constant X > J, we can rescale 
s — > sX/J and s ma x —> s max X/J, and with this rescaled 
s the requirement ||9 a /ii < J becomes satisfied. As- 
sume that H(0) has a known product ground state. Con- 
sider any observable O which is a product of operators 
supported on a single site. 

Our main result is that it is possible to compute the 
expectation value of O in the ground state of H(s max ) 
to any desired accuracy eo by an algorithm that takes a 
computational time T given by 



T 



exp 



[O(D°^)]O(N(J/AE)s max (N/e r), (1) 



on a classical computer, where the exponent to equals 

uj = o(\ti{D)J/Ae\ (2) 



For any operator P supported on a constant number of 
sites, n, we can write P as a sum over D 2n different prod- 
uct operators O. Therefore, the ability to approximate 
product operators implies the ability to approximate gen- 
eral operators on a constant number of sites. This result 
implies that physical quantities such as the ground state 
energy of H(s max ) can be approximated to within accu- 
racy l/poly(./V) in polynomial time on a classical com- 
puter. 

We rely heavily on the area law for one-dimensional 
systems with gapped Hamiltonians Q proven recently. 
Let tp°(s) denote the ground state of H(s). Given that 
H(s) satisfies the conditions above including J/AE be- 
ing 0(1), the area law implies that, for any e > 0, we can 



2 



approximate ip°(s) by a matrix product state ip m ps{s) 
such that 



\^ mps (s)-^<e 



(3) 



and such that tpmps(s) has bond dimension 

k = 0{N/ef exp[0(D ^)}. (4) 

The bond dimension is polynomial in N/e, but may be 
doubly exponentially large in J/ AE; this general upper 
bound applies to all one-dimensional systems while spe- 
cific cases in practice have not required such a large bond 
dimension. While the area law implies the existence of a 
good matrix product approximation to the ground state, 
it does not imply that we can efficiently find the cor- 
rect matrix product state. To solve this problem, we 
construct a sequence of matrix product approximations 
to the ground states along the entire path in parameter 
space, as we now explain (the explanation of how to do 
this given in Q was not correct, and we now give a full 
description of the correct construction). Given such a 
matrix product state, we can efficiently calculate expec- 
tation values of product operators such as O. 

The Algorithm — At s = 0, the ground state is as- 
sumed to be a known product state, and hence is a ma- 
trix product state of bond dimension k = 1. We break 
the adiabatic evolution for s = to s — s max into a se- 
quence of polynomially many discrete steps, such that s 
increases by a small amount in each step. Specifically, let 
a-max be the smallest integer larger than 4NJs max /AE, 
and let S = s max /a max so S < AE/ANJ. We break 
the adiabatic evolution into a max different discrete steps 
of size S, and we set s a = aS, for a — 0, 1, 2, a max . 
We have a max — 0(N (J / AE)s max ) . We fix a maximum 
bond dimension k maxi by setting k max = k for k given 
by Eq. (|4]) with e = l/poly(iV); the correct choice of e to 
obtain a given error is given after Eq. (f2"Tj) . Let ip® de- 
note the ground state of H(s a ). For a = 0, we represent 
the ground state ip® exactly as a matrix product state 
^™ ps . Our algorithm proceeds iteratively, taking a series 
of a max steps, such that on the a-th step it computes a 
matrix product state ip™ ps of bond dimension at most 
kmax which is a good approximation to ip®. We do this 
using the following algorithm: 

1: Initialize ^p™ ps to the known initial product state. 



2: For a = 1 to a r , 



do 



2a: Compute a matrix product state i/j^ from the 
state tpT-i as described in the section Improv- 
ing the Approximation. ijr a will be a good ap- 
proximation to "0a: but its bond dimension k' 
may be larger than k max by a factor polyno- 
mial in N and e. 

2b: Compute the state ?/'™ ps from the state ^ by 
a truncation procedure described in Trunca- 
tion Error Bounds. The state ipT ps wm be 



a matrix product state of bond dimension at 
most k max . 

We define the error after each step, e a , by 

l(VC ps >°)l 2 = i- ££ , (5) 

The correctness of the algorithm is based on the inductive 
assumption that after a — 1 steps, we have 



e a _i < min(AE/12NJ, 1/99), 



(G) 



so that is a good approximation to V'S-i- As we 

show in Eq. (fTT|) . since the difference \s a — s a -i| is less 
than or equal to AE/ANJ, the state ip™ P i is also a good 
approximation to -0°. We then use the result (fTT|) and 
the spectral gap and locality of the Hamiltonian to con- 
struct a state ip a which is a better approximation to ip® 
at the cost of increasing the bond dimension above k max 
as described two sections later. We then truncate ip a to 
a state ip™ ps with bond dimension k max , and we use the 
area law to bound errors in this truncation. 

Both the computational effort required to compute 
ifta and the bond dimension k' will be proportional to 
(N/e a )°( u \ Since it is possible, given a matrix product 
state of bond dimension k, to compute expectation values 
of observables which are products of single site operators 
in a time polynomial in fc, this shows the main result. 

Difference Between tpT-i an d Ya — Assume by induc- 
tion that after a — 1 steps we have a good approxima- 
tion ip™ P i to V'a-ii namely that Eq. ([6]) holds. Since 
I s a - Sq-i| < AE/4NJ, we have 



\H(s a ) - H{s a ^)\\ <AE/A. 



(7) 



From Eq. ©, 



(8) 

Let E® denote the ground state energy of H(s a ). We 
have \E° a - E° a _A < AE/A. Therefore, 

{^l-i,H{s a )ilP a _ 1 ) - E°< AE/2. (9) 

Since H(s a ) has a spectral gap AE, 

\(i>a^a-l)\ 2 > 1/2- (10) 

From Eq. ([5]), the angle Q\ between vectors i^T-i 
and Va-i obeys cos(6*i) 2 > 1 — e a _i. From Eq. (flQ)) . 
the angle # 2 between vectors ip^-i and obeys 
cos(# 2 ) 2 > 1/2. Therefore, we can bound the angle 
between vectors i/jT-i an< l ^2 % ^ - ^1 + ^2 with 
cos(0i + e 2 f = [cos(6>i)cos(6> 2 ) - sin(6»i)sin(6l 2 )] 2 > 
cos(6'i) 2 cos(6' 2 ) 2 -2cos(6'i)sin(6'i)cos(6l2)sin(6' 2 ) > (1 - 
e a _i)(l/2) — y/£a-i > 1/4. where the last inequal- 



ity follows because e a < 1/99 by (J( 

i(C- p r,^>i 2 = cos(0) 2 , 

KC- p r,^)i 2 >i/4 



Therefore, since 



(11) 
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Improving the Approximation — We now explain step 
2a, which uses a state 4>T-i which satisfies (fTTj) to 
construct, for any desired e, a state ip l a such that 
\(ipa, V'a)! 2 > 1 — e 7 with a bond dimensional and compu- 
tational cost of order (N/e)°( u \ In practice, we would 
prefer simpler constructions of the state tp^, as explained 
in the discussion. 

Let the states ip\ be a complete basis of eigenstates of 
H(s a ) with energies E%. From Eqs. (|7l6p and the fact 
that ||-ff(s a )|| < NJ, after a — 1 steps it is possible for 
the algorithm to compute E® to an accuracy of AE/A + 
(AE/12NJ)(NJ) = AE/3^. In the following, assume 
that the algorithm estimates E® to be zero (if this is not 
the case, shift the energies by the estimate of E®). 

Consider the state 

AE f 

<P a = -=j dtexp[-(AEt) 2 /2q]exp[iH(s a )t}^ 

= J2 H k a )^M-l{E k a /AEf/2] (^|C- P r) ) (12) 
fc 

where q is a number we choose below to be of order 
log(l/e), where e is the error estimate below in (JT3"jl . 
We will first show that the difference between (j) a /\<t>a\ 
and ipa is exponentially small in q and then we will show 
how to approximate cp a by a matrix product state, giv- 
ing the desired matrix product approximation to ip®. Let 
C-f = Aai& + B a ti where Vf) = 0. Note that by 
Eq. (HI]), we have \A a \ 2 > 1/4. Let cb a = A' a ^ a + B' a ^ 
with (ip®, <p^) =0. The idea of the integration over time 
is to approximately project onto the ground state; the 
projection of <p a on any state ip\ is reduced by a fac- 
tor of exp[— q(E^/AE) 2 /2] compared to the projection 
of VC-i ont0 iPa- Therefore, since \E°\ < AE/3, we 
have A' a > (1/4) exp(— q/18). However, since all other 
states have an energy at least AE above the ground state, 
and hence an energy at least 2AE /3 above zero, we have 
B' a < exp(— 2q/9). Thus, we can guarantee that the nor- 
malized state 4> a /\4>a\ is within error e' of ip^ for any e' 
by choosing q logarithmically large in N and e'. 

We now show how to approximate <p a by a matrix prod- 
uct state. First, replace the integral over all t between 
— oo and +oo by a finite integral, from t = —t max to 
t = +t max with t max = 99q/AE. The error in mak- 
ing this replacement is of order exp[— (AEt max ) 2 /2q] = 
exp(— 99 2 q/2). Also, replace the continuous integral 
over t by a discrete sum over different times ti. The 
bound on the operator norm of H(s a ), \\H(s a )\\ < NJ, 
implies that \dt exp[iH(s a )t]ip™ P i\ < NJ, and so we 
can approximate the integral by a sum with an er- 
ror e' using n sum = 0((t max AE/y / 2irq)(NJt max /e')) = 
0(y/qNJt max /e') terms in the sum. 

We now approximate the sum of these n sum states by 
a matrix product state with polynomial bond dimension. 
It was shown (llj using Lieb-Robinson boundsfH, [l]| 
that the state exp[iH(s a )ti]ip^ lP i can be approximated 



to error e' by a matrix product state with a bond di- 
mension exp(0(tiJ\nD))k max times a function of N/e' 
which grows slower than any power. Since ti is of or- 
der (1/AE) ln(N/e'), the state exp[iH(s a )ti}ip™v* can be 
approximated by a state with bond dimension of order 
(Nl e'^^hmn-r [ 1 4j] . The sum over n sum different matrix 
product states is still a matrix product state, with a bond 
dimension which is n sum times as large. Let the normal- 
ized sum of matrix product states be ip^. Therefore for 
any e we can choose q such that 



(13) 



and such that ip^ is a matrix product state with bond 
dimension polynomial in ./V and 1/e and such that com- 
puting ip l a requires only polynomial computational effort. 

Truncation Error Bounds — We now bound the trunca- 
tion error introduced in 2b. By the area law, the ground 
state ip° obeys 



\Xa 



(14) 



for some \a which is a matrix product state of bond 
dimension k, and thus from (|13I14|) we have 



|Xa-^| 2 <4e. 



(15) 



Pick any bond and do a Schmidt decomposition of ip l a 
across that bond, writing ip* = J2 a A(a)ipL{ct) <8> ^pn(a). 
Order the Schmidt coefficients so that |^4(a)| is decreas- 
ing as a increases. Then, from Eq. (fl5| . 



Y,\A(a)\ 2 <4e. 



(16) 



We now define ip™ ps - As shown in [15j, the bound on 
Schmidt coefficients in Eq. (fT!)]) implies that there exists 
a matrix product state ip™ ps with bond dimension k such 
that 



W-rtf <8(N-l)e 



(17) 



The construction in [15j is defined in terms of the ma- 
trices which are the parameters of the matrix product 
state; an equivalent construction is to define, for each 
i = 2,...,N, the operator Pi_N to project onto the k 
largest Schmidt coefficients of the reduced density ma- 
trix V'a on sites i, ...,N. Then, let ip™ ps be defined by 

C ps = z- 1 p N - hN p N - 2 , N ...p 2 ,Ni> t a , (18) 

By Eqs. PTT7D, 

|VC PS - ^a? < e + S(N - l)e + 2^/8(N - l)e. (19) 
Therefore, 

|(^,V2)| a >l-ea, (20) 
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with 



£a+l 



+ 8(N - l)e + 2y/8(N- l)e. 



(21) 



By picking e small enough, we can ensure that Q is true 
after a steps, given that it was true after a — 1 steps, and 
we can make t amax polynomially small at a polynomial 
computational cost. This completes the error estimates. 

Discussion — This work raises several natural ques- 
tions in the field of Hamiltonian complexity. First, what 
happens in more than one dimension? We do not know of 
a general proof of an area law in more than one dimen- 
sion, and it is quite conceivable that adiabatic simula- 
tion in two dimensions with a constant gap and a unique 
ground state is computationally universal. 

Second, the algorithm we have chosen is not very prac- 
tical. While the construction of the state ^* requires 
only polynomial time, and hence this work suffices to 
show in principle that adiabatic simulation with a con- 
stant gap is not useful in one dimension, in practice 
we prefer to construct matrix product states with as 
low a bond dimension as possible. This is why low or- 
der Trotter-Suzuki approaches are popular, such as the 
TEBD algorithm^!, [n}. Therefore, it is of great in- 
terest to prove that the following natural algorithm al- 
ways suffices to approximate the ground state: make 
the spacing s a — s a -i very small (polynomially small 
in N) so that the state if 1 ™-* is polynomially close to 
approximating the ground state tp®. Then, in step 2a, 
set tpQ = (1 — H(s a )/\\H(s a )\\)'4> 1 ^i, or perhaps instead 
tpl = exp(-i/ e . t , erl (s a )r)exp(-i? 0(irf (s a )T)V'™_ P i where r 
is some small quantity and -ff e i, e n,odd(s a ) represent the 
terms of H(s a ) on the even and odd bonds respectively. 
Finally, do the truncation step 2b as before. Such a proof 
would require a much more accurate analysis of the error 
in the truncation step. 

Finally, it is important to understand the role that 
the adiabatic evolution plays in our result. The area 
law implies the existence of a good matrix product state 
approximation to the ground state. This implies that 
the following decision problem is in NP: given a Hamil- 
tonian H with interaction strength J and local Hilbcrt 
space dimension both 0(1) and with nearest neighbor in- 
teractions, and given the promises that H has a unique 
ground state with inverse spectral gap which is 0(1), 
and that the ground state energy E is either < or 
> l/poly(7V), decide whether the ground state energy is 
indeed less than zero. Such a problem is in NP because 
the matrix product state guaranteed by the area law acts 
as a witness, in computer science language, or as a vari- 
ational state, in physics language. However, we have no 
guarantee that we can efficiently find such a state. In- 
deed, it has been shown that there are problems with a 
ground state which has a polynomial bond dimension as 
a matrix product state for which computing this matrix 



product state is NP-hard[18j. This result [18| holds for 
systems with an inverse polynomial spectral gap, rather 
than those with a constant gap as we consider here, so 
it is not certain what the difficulty is of finding matrix 
product states for the systems we consider, but it may 
indeed by a difficult problem. Thus, the ability to follow 
the state along the adiabatic change in the Hamiltonian is 
useful precisely because it allows one to be sure of locat- 
ing a good matrix product approximation to the ground 
state. 
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